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ABSTRACT 



We discuss optical associations, spectral energy distributions and photometric red- 
shifts for SWIRE sources in the ELAIS-N1 area and the Lockman Validation Field. The 
band-merged IRAC (3.6, 4.5, 5.8 and 8.0 fim) and MIPS (24, 70, 160 fim) data have 
been associated with optical UgriZ data from the INT Wide Field Survey in ELAIS-N1, 
and with our own optical Ugri data in Lockman- VF. Criteria for eliminating spurious 
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infrared sources and for carrying out star-quasar-galaxy separation are discussed, and 
statistics of the identification rate are given. 32 % of sources in the ELAIS-N1 field are 
found to be optically blank (to r = 23.5) and 16 % in Lockman-VF (to r = 25). 

The spectral energy distributions of selected ELAIS sources in Nl detected by 
SWIRE, most with spectroscopic redshifts, are modelled in terms of a simple set of 
galaxy and quasar templates in the optical and near infrared, and with a set of dust 
emission templates (cirrus, M82 starburst, Arp 220 starburst, and AGN dust torus) in 
the mid infrared. 

The optical data, together with the IRAC 3.6 and 4.5 fim data, have been used 
to determine photometric redshifts. For galaxies with known spectroscopic redshifts 
there is a notable improvement in the photometric redshift when the IRAC data are 
used, with a reduction in the rms scatter from 10 % in (1+z) to 7 %. While further 
spectroscopic data are needed to confirm this result, the prospect of determining good 
photometric redshifts for much of the SWIRE survey, expected to yield over 2 million 
extragalactic objects, is excellent. Some modifications to the optical templates were 
required in the previously uninvestigated wavelength region 2-5 /im. 

The photometric redshifts are used to derive the 3.6 and 24 /jm redshift distribution 
and to compare this with the predictions of models. For those sources with a clear 
mid infrared excess, relative to the galaxy starlight model used for the optical and near 
infrared, the mid and far infrared data are modelled in terms of the same dust emission 
templates (cirrus, M82, Arp 220, AGN dust torus). The proportions found of each 
template type are: cirrus 31 %, M82 29 %, Arp 220 10 %, AGN dust tori 29 %. The 
distribution of the different infrared sed types in the Li r /L opt versus L; Lr plane, where 
Li r and L opt are the infrared and optical bolometric luminosities, is discussed. 

There is an interesting population of luminous cool cirrus galaxies, with Lj r > L op t, 
implying a substantial dust optical depth. Galaxies with Arp220 like seds, of which 
there is a surprising preponderance compared with pre-existing source-count models, 
tend to have high ratios of infrared to optical bolometric luminosity, consistent with 
having very high extinction. There is also a high proportion of galaxies whose mid-ir 
seds are fitted by an AGN dust torus template (29%). Of these only 8 % of these are 
Type 1 AGN according to the optical-nir template fitting while 25% are fitted with 
galaxy templates in the optical-nir and have Lj r > L op t, so have to be Type 2 AGN. 
The remainder have Lj r < L opt so can be Seyferts, in which the optical AGN fails to be 
detected against the light of the host galaxy. The implied dust covering factor, > 75%, 
is much higher than that inferred for bright optically selected quasars. 

Subject headings: infrared: galaxies - galaxies: evolution - star formation - galaxies: 
starburst - cosmology: observations 
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1. Introduction 

Infrared wavelengths hold the key to understanding the evolution of galaxies. From the spec- 
trum of the extragalactic background radiation ( Hauser and Dwek 2001) we see that much of the 
light emitted by stars is absorbed by dust and reemitted at mid and far infrared wavelengths. Only 
by understanding infrared extragalactic populations can we hope to get a reliable census of the star 
formation history of galaxies and estimate the fraction of dust-obscured AGN. 

IRAS 12-25-60-100 /im colour-colour diagrams proved remarkably effective in identifying the 
main contributing components to the far infrared spectra of galaxies (Helou 1986, Rowan-Robinson 
and Crawford 1989). The key components were identified as infrared cirrus (emission from dust in 
the interstellar medium of the galaxy, illuminated by old stars and by quiescent star formation), 
starbursts (with the local starburst galaxy M82 as the prototype, generally associated with galaxy 
interactions), and AGN dust torus emission (predominantly at 3-30 //m in the rest-frame). Very 
high optical depth starbursts like Arp 220, probably associated with violent mergers, constituted 
a fourth component (Condon et al 1991, Rowan-Robinson and Efstathiou 1993). 

The Infrared Space Observatory broadly confirmed the reality of these infrared galaxy pop- 
ulations. Mid-infrared spectroscopy with ISO confirmed the distinct nature of these components 
and gave details of the mid-infrared PAH features (Xu et al 1998, Lutz et al 1998, Helou et al 
2000, Laurent et al 2000, Dale et al 2001, Lu et al 2003). The ELAIS survey (Rowan-Robinson et 
al 2004) provides the largest sample of galaxies with ISO photometry at 6.7, 15, 90 and 175 /im. 
Rowan-Robinson et al (2004) show optical and infrared spectral energy distributions of a sample 
of ELAIS galaxies, modelled in terms of a simple set of 4 templates, for each of which radiative 
transfer models exist. We have used the same four templates to fit the spectral energy distributions 
(seds) of SWIRE sources in this paper. 

The SPITZER SWIRE Legacy Survey is expected to detect over 2 million infrared galaxies 
in 49 sq deg of sky (Lonsdale et al 2003, 2004). This will be especially powerful for searches for 
rare objects and for systematic studies of the link between star formation history and large-scale 
stucture. To date we have examined IRAC 3-8 /im data and MIPS 24-160 /im data in a 0.3 sq 
deg validation area in our Lockman field (Lockman-VF), and in 6.5 sq deg of the ELAIS-N1 field. 
Although this is only a small sample of the SWIRE survey (~ 13%), we can already make a 
quantitative survey of the galaxy population. 

In the ELAIS-N1 area we can make a direct comparison with the ELAIS ISO survey (Rowan- 
Robinson et al 2004). Sources in common will have been observed in a total of 10 infrared bands 
between 3 and 200 /im, 15, 90 and 175 /im, from ISO and the 7 SPITZER bands. This allows us 
to carry out a detailed study of the spectral energy distributions (seds) of these galaxies. These 
studies validate our photometric redshift estimates and allow us to make estimates of the bolometric 
luminosity in different dust emission components. 

We present here a preliminary analysis of these seds and the characteristic luminosities of the 
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infrared galaxies in the two samples of sky. Polletta et al (2005 in prep) give a parallel discussion of 
the seds, redshifts and nature of SWIRE sources, with an emphasis on colour-colour diagram diag- 
nostics. Hatziminaoglou et al (2005) discuss SWIRE observations of SDSS quasars and Franceschini 
et al (2005) discuss SWIRE observations of Chandra X-ray sources in Nl. Surace et al (2005 in 
prep) discuss the SWIRE IRAC (3-8 /im) source-counts, Shupe et al (2005) discuss the SWIRE 24 
/im source-counts, and Oliver et al (2004) discuss the 3.6 ^m angular correlation function. 

The SPITZER telescope and mission are described by Werner et al (2004), the IRAC in- 
strument is described by Fazio et al (2004), and the MIPS instrument by Rieke et al (2004). 
Preliminary discussions of SPITZER galaxy populations have been given by Chary et al (2004) for 
the GOODS-SV survey, Eisenhardt et al (2004) for the IRAC Shallow Survey, Yan et al (2004) for 
the FLS Survey, and Lonsdale et al (2004) for the SWIRE Survey. The main thrust of these early 
discussions is the capability of SPITZER to detect star-forming galaxies at z > 1. However these 
discussions were hampered by the absence of redshift information. In this paper we use photometric 
redshifts to give a fuller picture of the infrared galaxy populations. 

A model with f2o = 0.3, Ao = 0.7 and a Hubble constant of 72 km/s/Mpc is used throughout. 

2. Optical associations 

Our starting point is the band-merged SWIRE 3.6-24^m point-source catalogues (Surace et 
al 2004) , which consist of sources detected with an SNR of at least 5 in one or more IRAC bands 
and their 24 /im associations with an SNR of at least 3. The astrometric accuracy of the IRAC 
sources appears to be very good, better than ±0.5 arcsec (2-<r). Aperture fluxes are used for 
24 jj,m and fainter IRAC sources, but Kron fluxes, which give a better estimate of the flux from 
extended sources, have been used for the IRAC bands if S(3.Q)Kron > 1 mJy. Our goal is to 
define samples of good-quality SWIRE sources with reliable optical identifications and sufficient 
photometric bands for photometric redshift estimation. We search for optical associations of the 
SWIRE sources, using the band-merged WFS survey UgriZ data for the Nl area (Babbedge et 
al 2005, in preparation), complete to r ~ 23.5, and our Ugri survey data for Lockman-VF (Siana 
et al 2005, in preparation), complete to r ~ 25, using a search radius of 1.5". Analysis of the 
identification rate as a function of search radius shows that an increasing fraction of associations 
at radii larger than this are spurious (Surace et al 2004). Where a SWIRE source picks up more 
than one optical association, the association with the most optical bands is selected first (to select 
against a number of spurious single band optical sources, and to favour brighter sources) and then 
the nearest if more than one have the same number of optical bands. Where two infrared sources 
closer than 0.5 arcsec acquire the same optical association we have kept only the nearer, since 
we suspect the majority of these close pairs may be processing artefacts. However there are some 
genuinely blended sources which will require a more sophisticated treatment. After this elimination 
of duplicates there are 16149 band-merged (IRAC+MIPS24) sources in the 0.3 sq deg Lockman 
area with optical associations in at least one optical band, and 204866 band-merged sources in the 
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6.5 sq deg of ELAIS-N1. Examination of band-merged sources with no optical association showed 
that some of these are spurious, generated for example by bright stars (see discussion by Surace et 
al 2004). The requirement that appears to eliminate spurious sources effectively and give a high 
reliability catalogue is that the source be detected at both 3.6 and 4.5 fim at SNR > 5 (Surace 
et al 2004). This yields 96155 blank fields in Nl and 2991 in Lockman-VF, so that the identified 
sources represent 68 and 84 % of the total. Restricting to sources detected at 24 /im (38601 in 
ELAIS-N1), 35 % are optically blank (to r = 23.5). For comparison Yan et al (2004) found 38 % of 
24 ^m sources in the FLS Survey have r > 23.0 and 17 % have r > 25.5. SWIRE blank fields will 
be discussed further in a subsequent paper (Lonsdale et al 2004, in preparation). Table 2 gives a 
summary of the statistics of the SWIRE catalogue in the ELAIS-N1 area. 

Finding charts for ELAIS optical identifications and a summary of optical properties can be 
found in Gonzalez-Solares et al (2004) and at http://www.ast.cam.ac.uk/ eglez/eid/ 

3. Star-galaxy separation 

Optical objects classified as extended by Sextractor are assumed to be galaxies. Those classified 
as 'stellar' may be stars, quasars or distant galaxies. As an initial cut at removing the stars, we 
removed objects with r < 23 classified as stellar in r. However some of these will be quasars, so 
we defined an initial 'quasar' subset as (1) objects with an 8 /im/r-band colour cooler than a 2000 
K blackbody (due to presence of an AGN dust torus), or (2) objects whose UgriZ 3.6, 4.5 /im seds 
are best fitted with an AGN rather than a galaxy template. 

Figures 1 and 2 show the 3.6, r, i colour-colour diagrams for the 'stars' and 'quasars' defined 
in this way, for tile 3-2 in Nl (the SWIRE Nl data was analyzed in 0.8 sq deg tiles, and for some 
figures in this paper it is clearer to show results only for a single tile). Clearly the 'stars' still 
include a significant number of extragalactic objects, which we can recover by keeping those above 
the straight line locus shown. Similarly the 'quasars' still have a stellar contamination which can 
be removed by omitting objects below the line shown. To show how effective this procedure is, 
Figure 3 shows the g-r-i colour diagram for our final stellar catalogue in Nl, with almost all objects 
following the stellar colour locus. We find that 13 % of Nl sources and 6 % of Lockman-VF sources 
are stars. A similar star-galaxy separation criterion has been proposed by Eisenhardt et al (2004). 

A few bright stars do not satisfy the stellar rejection criterion because of saturation in the 
optical bands. They are easily recognized as having r < 16, stellar or saturated flags, photometric 
redshifts close to zero and very poor x 2 , and have been reclassified 




Fig. 1. — 3.6 /xm-r-i colour-colour diagram for preliminary 'stars' catalogue in Nl (tile 3-2). The 
line has equation 3.2(r-i+0.5) = 3.5 (log S(3.6)/S(r) +1). 
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Fig. 2. — 3.6 /xm-r-i colour-colour diagram for preliminary 'quasars' catalogue in Nl (tile 3-2). 



Fig. 3. — g-r-i colour-colour diagram for final stars catalogue in Nl, with 17 < r < 22. Fig 3 
available at http://astro.ic.ac.uk/ mrr/ 
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4. Spectral energy distributions 

To model the seds of SWIRE sources we have used an approach similar to that of Rowan- 
Robinson et al (2004) for ELAIS sources. The optical data and near infrared data to 4.5 fim are 
fitted with one of the 8 optical galaxy templates used in the photometric redshift code of Rowan- 
Robinson et al (2003), 6 galaxy templates (E, Sab, Sbc, Scd, Sdm, sb) and 2 AGN templates 
(Rowan- Robinson et al 2004) . The mid and far infrared data are fitted by a mixture of the 4 infrared 
templates used by Rowan-Robinson (2001) in models for infrared and submillimetre source-counts: 
cirrus, M82 starburst, Arp 220 starburst and AGN dust torus. For each of these templates we 
have detailed radiative transfer models (cirrus: Eftstathiou and Rowan-Robinson 2003, M82 and 
Arp 220 starbursts: Eftstathiou et al 2000, AGN dust tori: Rowan-Robinson 1995, Efstathiou and 
Rowan-Robinson 1995). The approach here is to try to understand the overall SPITZER galaxy 
population. Better fits to individual galaxies could be obtained by allowing variation of the many 
parameters in the radiative transfer codes from which these templates are taken and this will be 
explored in subsequent work. 

With the first tranche of SWIRE data in the ELAIS-N1 region, we are able to fill in a major 
wavelength gap in these seds and show for the first time high quality infrared seds from 0.36-175 
jum, providing a strong test of these models. In figures 4-6 we show sed fits for 30 sources with 
spectroscopic redshifts which are detected in most of the ISO and SWIRE bands (including MIPS 
70 and 160 /im). The ELAIS Final Catalogue was associated with the SWIRE 3.6-24 fim catalogue 
using a search radius of 2.5 arcsec. The MIPS 70 and 160 /im fluxes were extracted from mosaiced 
maps by Donovan Domingue using the Starfinder code (Diolaiti et al 2000). 70 and 160 p source- 
lists were associated with the SWIRE sources with 24 ^m detections using a search radius of 7 
arcsec at 70 fim and 14 arcsec at 160 fim. Histograms of source separation show an increasing 
fraction of spurious associations at radii greater than these. Where 70 or 160 fi sources acquire 
more than one 24 /jm association, the flux was always assigned to the brighter 24 fim source. In 
some cases, especially at 160 ^m, the flux should probably be assigned to more than one shorter 
wavelength source, so the 160 fim flux will be overestimated. For a few sources which are extended 
at 70 or 160 fim, these fluxes may be underestimated. 916 sources were associated in this way with 
70 or 160 fim detections, with about 35 % of these being ELAIS sources, within the region of overlap 
of the two surveys. The numbers of sources in the fully bandmerged 3.6-175 /im SWIRE+ELAIS 
catalogue with 10, 9, 8 or 7 infrared bands out of 3.6, 4.5, 5.8, 8.0, 15.0, 24.0 70, 90, 160, 175 ^m 
is 14, 44, 93, 332, respectively. 

Figure 4 shows seds for 10 galaxies with spectroscopic redshifts and which are detected at 
UgriZ, JHK (primarily from 2MASS), 3.6, 4.5, 5.8 8.0, 15, 24, 70, 90, 160 and 175 /im, and which 
were characterized by Rowan-Robinson et al (2004) as cirrus galaxies, ie normal spirals with mid 
and far infrared emission from interstellar dust illuminated by the general stellar radiation field. 

Table 1 gives the parameters of the model fits: SPITZER RA and dec, ELAIS name, redshift 
(bracketed if photometric), n se( i (infrared sed type: 1 = cirrus, 2 = M82 starburst, 3 = Arp220 
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Fig. 4. — Spectral energy distributions for 10 galaxies with spectroscopic redshift from the ELAIS 
survey, characterized by Rowan-Robinson et al (2004) as cirrus galaxies, with 10 infrared spectral 
bands at A > 3/«n. The SPITZER data fit well onto the cirrus template model of Rowan-Robinson 
(2001, see also Efstathiou and Rowan-Robinson 2003), and there is good agreement between the 
ISO and SPITZER fluxes. Dotted curves show the optical-nir and cirrus templates, solid curve the 
total predicted sed. Data are the INT Wide Field Survey UgriZ, JHK (generally 2MASS), ISO 
15,90, 175 fim, IRAC 3.6, 4.5, 5.8, 8.0 ftm, and MIPS 24, 70, 160 fim. Error bars (typically ±0.04 
dex, corresponding to a 10% photometric uncertainty) are smaller than the size of the symbols. 
Details of the sources in Figs 4-6 are given in Table 1, sequentially starting from the bottom of Fig 
4. 
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Fig. 5. — Spectral energy distributions for 10 galaxies, 9 with spectroscopic redshifts, from the 
ELAIS survey. 3 are luminous cirrus galaxies, while 7 are fitted by a combination of an M82 
starburst and cirrus. The optical-nir and M82 templates are shown as dotted curves, the cirrus 
template as a broken curve, and the total predicted sed as a solid curve. 
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Fig. 6. — Spectral energy distributions for 5 galaxies from the ELAIS survey characterised as Arp 
220 starbursts, 3 with spectroscopic redshifts, and 5 quasars with spectroscopic redshifts, with 
infrared seds modelled with an AGN dust torus and, in two cases, an additional M82 starburst 
component. The optical-nir and Arp 220 templates are shown as dotted curves, the AGN dust 
torus as a long-dashed curve, the M82 starburst as a short-dashed curve, and the total predicted 
sed as a solid curve. 
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starburst, 4= AGN dust torus), log 10 bolometric infrared luminosity Lj r , the vertical offset constant 
C used in plotting the seds, the optical-nir sed type ntyp (1 = E, 2 = Sab, 3 = Sbc, 4 = Scd, 5 = 
Sdm, 6 = sb, 7 and 8 = QSOs, see Rowan-Robinson 2003), Ay, logio bolometric optical luminosity 
L opt . The data at wavelengths < 5 /im are fitted with the optical galaxy templates used by Rowan- 
Robinson et al (2004) in their photometric redshift estimation (see section 5). The spectroscopic 
redshifts are in the range 0.03-0.27. The cirrus template used at wavelengths > 5 fim is taken 
from the sequence of radiative transfer models by Efstathiou and Rowan- Robinson (2003). The 
SPITZER data fit remarkably well onto the models determined solely from ISO, optical data and 
JHK data, demonstrating the consistency of the IRAC and MIPS calibrations. It is remarkable how 
well a single cirrus template fits these relatively quiescent galaxies (Lj r < 2.10 11 L Q ). A few galaxies 
(not shown here) will require a cooler cirrus template. There is a reasonable consistency between 
the ISO 90 and 175 fim data and the MIPS 70 and 160 ^m data. For example for cirrus sources in 
common it is found that log 10 (S (90) / S (70)) = 0.27 + 0.17 and log w (S '(175) /S '(160)) = 0.02 ±0.12, 
so the consistency of both ISO and SPITZER calibrations is supported. 

Figure 5 shows seds for 10 ELAIS-N1 galaxies which are characterized as higher-luminosity 
cirrus galaxies or M82 starbursts, all except one detected in at least two of the bands 15, 90 and 175 
fim. The M82 template is taken from the sequence of radiative transfer models by Efstathiou et al 
(2000). The galaxies have spectroscopic redshifts in the range 0.2-0.9 . The three 'cirrus' galaxies 
are very luminous, with Lj r > 3.10 n L Q and have Lj r significantly larger than L opt , implying a 
substantial optical depth in dust. One (241.67583+55.83316) is ultraluminous (L ir = 2.4.1O 12 L ), 
based on a photometric redshift of 0.479. Mid infrared spectroscopy and submillimetre photometry 
are needed to confirm the existence of this class of cool luminous galaxies, whose existence is 
predicted by the counts model of Rowan-Robinson (2001) and for which direct evidence is presented 
by Efstathiou and Rowan- Robinson (2003), Chapman et al (2003), Kaviani et al (2003) and Rowan- 
Robinson et al (2004). These are discussed further in section 7. 

Figure 6 shows seds for 10 ELAIS-N1 galaxies, comprising 5 Arp 220 (higher optical depth) 
starbursts, and 5 AGN dust tori. Two of the latter also show evidence of a starburst component at 
longer wavelengths. The Arp 220 template is taken from the sequence of radiative transfer models 
by Efstathiou et al (2000) and the AGN dust torus template is taken from the radiative transfer 
model of Rowan-Robinson (1995). The objects have spectroscopic redshifts in the range 0.06-3.6 
(two of the Arp 220 objects have photometric redshifts only). The fits for 160901+541808 and 
161005.8+541029 are worse than for other objects plotted at A > 60/xm. 

The striking features that emerge from this modelling are (1) that the seds can at least broadly 
be understood in terms of a small number of infrared templates, (2) the emergence of the class 
of luminous cool galaxies with Lj r > L opt , (3) the significant incidence of Arp 220 type seds in 
the survey, which was not predicted by source count models (eg Rowan- Robinson 2001). Radiative 
transfer codes have several parameters and even better fits could be found by using the full range of 
models of Efstathiou et al (2000), Eftstathiou and Rowan-Robinson (2003). This will be explored 
by Efstathiou et al (2004, in preparation). 
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In fitting these seds, and in the photometric redshift analysis described in the next section, it 
became apparent that the predictions of the optical-nir templates used by Rowan-Robinson (2003) 
at wavelengths 2-5 /um, where they had not previously been tested, were too low. Figures 7 and 8 
show a comparison of data for galaxies with spectroscopic redshifts and which are best fitted by E 
and Scd templates, respectively. Stellar synthesis fits to these seds (eg Babbedge et al 2004) will 
require a larger population of old low-mass stars. 

5. Photometric redshifts in ELAIS-N1 and Lockman-VF 

We have been able to carry out a photometric redshift analysis for the two SWIRE samples, 
the 0.3 sq deg area of Lockman-VF, in which we have moderately deep Ugri optical data (to r 
~ 25), and the 6.5 sq deg of the Nl area, in which we have the re-reduced WFS UgriZ data 
(to r ~ 23.5). Full analyses of the photometric data in these fields will be reported by Siana et 
al (2005, in preparation) and Babbedge et al (2005, in preparation). Figure 9 shows a plot of 
logio(S(3.Q)) versus r for the Lockman-VF field, which illustrates that the optical selection starts 
to produce serious incompleteness below S(3.6) ~ 40 fiJj for an optical survey complete to r = 
23.5 (as in ELAIS-N1), and below S(3.6) ~ 10 [iJy for an optical survey complete to r = 25 (as in 
the Lockman-VF). Figure 10 shows logio(S(24)) versus r for a 0.8 sq degree sub-area of Nl, with 
symbols and colours denoting the different infrared template types, discussed in section 7. Optical 
incompleteness becomes serious for 24 /im sources below 200 /uJy for a survey complete to r = 23.5. 

In this paper we have analyzed the catalogues of SWIRE sources with optical associations with 
the code of Rowan-Robinson (2003) which uses 6 optical galaxy templates and 2 AGN templates (as 
described in Rowan- Robinson et al 2004). The galaxy templates were based on empirical seds by 
Yoshii and Takahara (1988) and Calzetti and Kinney (1992) but have been subsequently modelled 
with a stellar synthesis code by Babbedge et al (2004). Parallel work is being carried out with 
IMPz, a more detailed version of this code described by Babbedge et al 2004 (Babbedge et al 2005, 
in preparation), and with a modified version of the hyper-z code (Bolzonella et al 2000) developed 
by M.Polletta (Polletta et al 2005 in preparation). Agreement between these three approaches 
applied to ELAIS data was found to be good by Babbedge et al (2004). 

We set a further condition on the optically associated SWIRE sources for entry to the photo- 
metric redshift code, namely that the source be detected at 3.6 /um, in the r-band, and in at least 
two other bands. Rowan-Robinson (2003) found that four bands was the minimum for effective 
photometric redshift estimation. This reduced the sample of optically identified SWIRE sources 
(after elimination of stars) to 133645 in Nl and 12664 in Lockman-VF. For only 6 % of these 4-band 
optically associated extragalactic SWIRE sources in Nl does the code fail to produce a redshift, 
with the reduced \ 2 > 10, and 12 % in Lockman-VF. To quantify the goodness of fit, figure 14 
shows the histogram of reduced \ 2 f° r Nl galaxies, which is consistent with the templates being 
used providing a reasonable fit to the optical-near-infrared data. 
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Fig. 7. — Comparison of normalized photometric data for SWIRE-ELAIS galaxies with spectro- 
scopic redshifts whose seds are classified as E with optical-nir sed template for elliptical galaxies 
used in the photometric redshift code, illustrating the need to refine the templates at wavelengths 
> 2fim. 
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Fig. 8. — Comparison of normalized photometric data for SWIRE-ELAIS galaxies with spectro- 
scopic redshifts whose seds are classified as Scd with optical-nir sed template for Scd galaxies used in 
the photometric redshift code, illustrating the need to refine the templates at wavelengths > 2/xm. 
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Fig. 9. — logioS(3.6) versus r for identified SWIRE sources in Lockman-VF, with different symbols 
for infrared template types: filled circles (black): cirrus, filled triangles (red): M82 starburst, filled 
squares (green): Arp220, crosses (blue): AGN dust tori. Sources with no infrared excess or excess in 
only a single band are shown as yellow open circles. Vertical lines indicate the optical completeness 
of the WFS and Lockman-VF surveys. 
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Fig. 10. — logioS(24:) versus r for identified SWIRE sources in Nl (tile 2-2), with different symbols 
for infrared template types: filled circles (black): cirrus, filled triangles (red): M82 starburst, open 
squares (green): Arp220, crosses (blue): AGN dust tori. 
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For a sample of the brighter Lockman-VF galaxies (the WIYN survey, Owen et al 2004, in 
preparation), and for 25% of the galaxies in Nl which are common to the ELAIS survey (Perez- 
Fournon et al 2004, Serjeant et al 2004, summarized in Rowan- Robinson et al 2004), we have 
spectroscopic redshifts and can make a direct comparison with the derived photometric redshifts. 
An important innovation offered by the SPITZER data is that we can use the 3.6 and 4.5 //m fluxes 
in the photometric redshift fits. Dust emission in cirrus and starburst galaxies generally starts to be 
a significant contributor to the total flux at around 5 /im in the rest frame. We find that inclusion 
of extinction, characterized by Ay, is essential for a sensible redshift solution for SWIRE galaxies, 
some of which appear to be dusty, with Ay ranging up to 2. We also find it necessary to allow 
modest extinction for the QSOs. We restricted Ay for QSOs to < 0.3, since allowing a greater 
extinction range lets in many more aliases for normal galaxies. This may result in some heavily 
reddened quasars not achieving a correct fit. Figure 11 shows a comparison of spectroscopic and 
photometric redshifts for 120 galaxies in the SWIRE-N1 area, derived from UgriZJHK data plus 
SWIRE 3.6 and 4.5 /im fluxes. The rms deviation of logio(l + z) for galaxies is 6.9 %, a significant 
improvement. For comparison typical rms values found by Rowan- Robinson (2003) and Rowan- 
Robinson et al (2004) from UgriZ, JHK data alone were 10 %. Figure 12 shows a comparison 
of spectroscopic and photometric redshifts for 289 galaxies in the WIYN survey in Lockman-VF 
derived from Ugri photometry, plus 3.6 and 4.5 ^m fluxes. The rms deviation of log\o{l + z) is 8.4 
%, again a significant improvement over determinations from Ugri data alone (the corresponding 
plot when 3.6 and 4.5 /^m data are not used is shown in Fig 13). Use of IRAC data in photometric 
redshift determination has also reduced the number of outliers. Redshift estimates for AGN can 
be affected by optical variability since optical photometry for different bands may have been taken 
at different epochs (Afonso-Luis et al 2004). A fuller discussion of photometric redshifts of SWIRE 
galaxies will be given by Babbedge et al (2005, in preparation). 

The distribution of optical-ir sed types for the 126193 galaxies in Nl for which we have pho- 
tometric redshift estimates is given in Table 2. The number of ellipticals my be inflated by objects 
similar to Arp 220, whose optical sed looks very like an elliptical because of the heavy obscuration 
of the newly- forming stars. 

6. Redshift distributions 

Figure 15 shows the redshift distributions derived in this way for SWIRE-N1 sources with 
S(3.6) > 40 jtxJy, above which flux the optical identifications (to r ~ 23.5) are, according to Fig 
9, relatively complete, with a breakdown into elliptical, spiral + starburst and quasar seds based 
on the photometric redshift fits. For comparison we show in the top panel the predictions of the 
Rowan-Robinson (2001), Xu et al (2003) and Pozzi et al (2004) models. Ellipticals cut off sharply 
at around z ~ 1.0, and spirals at z ~ 1.3 

Figure 16 shows the redshift distribution for Lockman-VF sources with S(3.6) > 10 /iJy, for 
which the optical identifications (to r ~ 25) are relatively complete. In the top panel we show 



0.2 0.4 0.6 0.8 

lgio(l + z sp ) 

Fig. 11. — Photometric redshifts in ELAIS-N1 using UgriZ, JHK, and SWIRE 3.6 and 4.5 
data, versus spectroscopic redshift. Filled circles: galaxies, filled triangles: QSOs. The ±2cr loci 
are indicated. Spectroscopic redshifts are from Perez-Fournon et al (2004, in preparation), Serjeant 
et (2004) (see also Rowan- Robinson et al 2004), Hatziminaoglou et al 2004. 
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Fig. 12. — Photometric redshifts for the WIYN survey using Ugri and SWIRE 3.6 and 4.5 fim 
data, versus spectroscopic redshift. Filled circles: galaxies, filled triangles: QSOs. The ±2a loci 
are indicated. 
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Fig. 13. — Photometric redshifts for the WIYN survey using Ugri data only, versus spectroscopic 
redshift. Filled circles: galaxies, filled triangles: QSOs. 
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Fig. 14. — Histogram of reduced x 2 for the photometric redshift fits. 
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predictions for the same 3 models for 5(3.6) > lOfiJy. Again the ellipticals cut off at z ~ 1.4, but 
a tail of star- forming galaxies is seen to z ~ 3. 

The photometric estimates of redshift for AGN are more uncertain than those for galaxies, 
due to aliassing problems and to the effect of variability on photometry, but the code is effective 
at identifying Type 1 AGN from the optical and near ir data. For some quasars there is significant 
torus dust emission in the 3.6 and 4.5 //m bands, so inclusion of these bands in photometric redshift 
determination can make the fit worse rather than better. We have therefore omitted the 3.6 and 
4.5 fim bands if S(3.6)/S(r) > 4. Note that only 5 % of SWIRE sources are identified by the 
photometric redshift code as Type 1 AGN. 

Extrapolating from Figs 9 and 10, we estimate that no more than 5-10 % of sources are missing 
from Figs 15 and 16 due to the optical incompleteness of these two surveys. However the actual 
proportions of sources which are either optically blank or have too few optical bands for photometric 
redshift determination is much higher than this, 25-35 %. Further investigation is needed to confirm 
the reality of these sources and their nature. Potentially they could substantially shift the shape of 
the redshift distribution towards higher redshifts. Further discussion of blank field sources is given 
in section 7. 



7. Bolometric infrared and optical luminosities 

For sources detected at 70 or 160 /im, or with an infrared excess at 4.5-24 /xm, relative to the 
template used for photometric redshift fitting, in at least two bands (one of which we require to be 
8 or 24 jum), we have determined the best-fitting out of cirrus, M82 starburst, Arp220 starburst 
or AGN dust torus infrared templates (as defined in section 3). From the experience of modeling 
IRAS sources (eg Rowan- Robinson and Crawford 1989), confirmed by the sed modeling described 
in section 3, we allow (a) a mixture of cirrus and M82 templates, (b) a mixture of AGN dust 
torus and M82 templates, or (c) an Arp 220 template (using a mixture of AGN dust torus plus 
A220 template did not improve the x 2 distribution). We classify all galaxies with a parameter 
n se d which is 1 for for an sed dominated by cirrus, 2 for sed dominated by M82 starburst, 3 for 
Arp 220 starburst, 4 for sed dominated by an AGN dust torus, 5 if there is an infrared excess in 
only a single band (in this case we can not say what the appropriate infrared template is), and 6 
for no infrared excess (bulge or elliptical dominated, or simply too faint for dust emission to be 
detected). Several figures in this paper have been colour-coded using this parameter. Table 2 gives 
the break down of the Nl sources with photometric redshift by their mid-ir (4.5-24 fj,m) template 
type. The proportion of cirrus galaxies (31%) agrees well with the figure given by Yan et al (2004) 
for the Spitzer First Look Survey (30%). Note the high proportion of galaxies whose mid-ir seds 
are fitted by an AGN dust torus template (29%). Only 8 % of these are Type 1 AGN according to 
the optical-nir template fitting. 25% are fitted with galaxy templates in the optical-nir and have 
Li r > L opt and so have to be Type 2 AGN. The remainder have Lj r < L opt so can be Seyferts, 
in which the optical AGN fails to be detected against the light of the host galaxy. The implied 
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Fig. 15. — Photometric redshift histogram for SWIRE ELAIS-N1 sources with good 4-band optical 
IDs, and S(3.6) > 40 fj,Jy. Top panel: all sources, dotted curve: prediction of Rowan-Robinson 
(2001), short-dashed curve: prediction of Pozzi et al (2004), long-dashed curve: prediction of 
Xu et al (2003, model E); lower panels: separate breakdown of contributions of ellipticals, spirals 
(Sab+Sbc+Scd), starbursts (Sdm+sb), and quasars. The histogram for quasars has been multiplied 
by 20. 
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Fig. 16. — Redshift histogram for SWIRE Lockman-VF sources with good 4-band optical IDs, and 
S(3.6) > 10 fJ,Jy. Top panel: all sources, dotted curve: prediction of Rowan- Robinson (2001), short- 
dashed curve: prediction of Pozzi et al (2004), long-dashed curve: prediction of Xu et al (2003, 
model E); lower panels: separate breakdown of contributions of ellipticals, spirals (Sab+Sbc+Scd), 
starbursts (Sdm+sb), and quasars. The histogram for quasars has been multiplied by 10. 
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dust covering factor, > 75%, is much higher than that inferred for bright optically selected quasars, 
~ 30% (Rowan-Robinson 1995). 

We can estimate the bolometric luminosity corresponding to the infrared template and to the 
optical template used for photometric redshift determination. Figure 17 shows the ratio of bolomet- 
ric infrared to optical luminosity, lgio(Li r /L op t), versus bolometric infrared luminosity, lg±oLi r for 
916 galaxies detected in Nl at either 70 or 160 /im. (or both), with different symbols/colours for the 
different infrared sed types. For star-forming galaxies the parameter Li r /L opt can be interpreted 
as approximately 10~ 9 M^/M^(yr~ 1 , since Lj r ~ 10 10 M* and L opt ~ 10M*, i.e. as 10~ 9 r~ 1 , where 
r is the time-scale in yrs to accumulate the present stellar mass, forming stars at the current rate 
(M^/M* = t~ 1 is the specific star formation rate). For galaxies with L{ r < 1O 1O L , a significant 
fraction of the infrared emission is reemission of starlight absorbed by (optically thin) interstellar 
dust, so Li r /L pt should be interpreted as the optical depth of the interstellar dust. Many very 
low values of Li r /L opt (< 0.2) are due to elliptical galaxies with a small amount of star-formation. 
Fig 18 shows the same plot for 2913 galaxies based on 4.5-24 fim data from a 0.8 sq deg sub-area 
of Nl (tile 2-2). These figures contains a wealth of information. The cirrus galaxies at low in- 
frared luminosities (< 10 11 Lq) typically have values of lgio(Li r / L op t) of -0.2, in agreement with 
the values found for nearby IRAS galaxies (Rowan- Robinson et al 1986). However there is the 
interesting population of luminous and ultraluminous cirrus galaxies, with L,- ir > L opt , discussed in 
section 3. 7 % of the cirrus galaxies identified in the Nl area have Lj r > 3.10 11 L Q and 2 % have 
Li r > 10 12 L Q . The implications are that (1) the quiescent phase of star formation was significantly 
more luminous in the past (as assumed in the count models of Rowan-Robinson 2001), (2) the 
dust opacity of the interstellar medium in galaxies was higher at z ~ 1, as expected from galaxy 
models with star-formation histories that peak at z = 1-2 (Pei et al 1999, Calzetta and Heckaman 
1999, Rowan- Robinson 2003). The galaxies with Arp 220 templates tend to have high values of 
Lir /L opt , consistent with the idea that they are very dusty. M82 type starbursts are less extreme 
and range over rather similar values of L{ r and Li r /L opt to the cirrus galaxies. Figure 19 shows a 
histogram of reduced x 2 f° r the infrared sed template fitting, for an assumed typical error of 0.06 
dex (corresponding to ±15%), which gives a measure of the goodness of fit of the ir templates. 

For galaxies detected at 70 or 160 /jm we can compare the infrared template fits and bolometric 
luminosity estimates derived from 4.5-160 //m data, and from 4.5-24 fim only. For cirrus galaxies 
the agreement is good but for other types it is not so good and for all sources the rms value 
Z<?io(-k&oZ,4.5-i6o/-k&oJ,4.5-24) is 0.30, corresponding to a factor 2 uncertainty in L^i estimated from 
4.5-24 fim only. This is unfortunate because we have 70 or 160 /j,m data for only a small percentage 
of the SWIRE galaxies with photometric redshifts. 

To summarise (see Table 2), for 63 % of the 126193 SWIRE Nl galaxies with photometric 
redshifts, we have no information at all on the dust emission (no significant infrared excess relative 
to the starlight/QSO models), For 36 % we have an approximate estimate based on 4.5-24 /im data 
which gives the bolometric luminosity only to an accuracy of a factor of 2, and for just 0.7 % of 
galaxies do we have a more accurate estimate based on longer wavelength data. 
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Fig. 17. — Ratio of bolometric infrared to optical luminosity, lgio(Li r /L op t), versus bolometric 
infrared luminosity, lgioLi r , for all sources in Nl detected at either 70 or 160 fim. Filled circles 
(black): galaxies with cirrus seds, filled triangles (red): M82 starbursts, open squares (green): 
A220 starbursts, crosses (blue): AGN dust tori. Right-hand scale gives values of the specific star- 
formation rate, lgio(M*/M*), (yrs~ 1 ). 
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Fig. 18. — Ratio of bolometric infrared to optical luminosity, lgio(Li r /L op t), versus bolometric 
infrared luminosity, IgioLfr, for SWIRE-N1 (tile 2-2) sources with photometric redshifts and at least 
2 bands at A > 4.5 fim with excess infrared emission above the optical galaxy template. Filled circles 
(black): galaxies with cirrus seds, filled triangles (red): M82 starbursts, open squares (green): A220 
starbursts, crosses (blue): AGN dust tori. Right-hand scale gives values of Z#io(-M*/-M*), (yrs^ 1 ). 
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Fig. 19. — Histogram of reduced \ 2 f° r the ir template fits, assuming typical error of 0.06 dex. 
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The code of Rowan-Robinson (2003) imposes an upper limit on the B band luminosity for 
starlight, which translates approximately to a limit on the optical bolometric luminosity (and 
stellar mass), which manifests itself as the sloping cutoff to the right of Fig 17. To the right of this 
limit can be seen the Type 1 quasars, some of which are best fitted by AGN dust tori, but others 
prefer M82 starbursts. While type 1 AGN would be expected to lie below the L ir = L opt line, since 
only part of the optical light we see is being intercepted by the dust torus, many of the galaxies 
whose mid and far infrared seds are dominated by AGN dust tori have Lj r > L op t, consistent with 
being Type 2 AGN. Further discussion of SWIRE quasars is given by Hatziminaoglou et al (2004). 

Quite a few objects are seen to have extremely high infrared luminosities, in the hyper luminous 
class (> 10 13 L Q ), most of them associated with quasars at high redshifts. Clearly spectroscopy is 
needed to determine these redshifts more accurately. 

Since the sources plotted in Fig 18 are mostly detected at 24 /mi it is of interest to consider 
whether PAH features are affecting the detectability of 24 /mi sources. Figure 20 shows redshift 
distributions for 24 /mi broken down by infrared sed type. Cirrus sources, and M82 starburst show 
a strong cutoff at z ~ 1.4, which may be partly due to the onset of the 10 /im silicate absorption 
feature. For Arp 220 sources, the cutoff starts at lower redshift (~ 1), consistent with their deeper 
and broader silicate absorption. A secondary peak for M82 starbursts (and possibly cirrus sources) 
at z = 0.8-1.4 can be interpreted as sources having their detectability at 24 /mi enhanced by the 
10-12 /im PAH feature. Sources with an infrared excess in a single band have a broad distribution 
between z = 0.2 - 1.4, and are consistent with being M82 starbursts. The optically blank 24 /mi 
sources would be likely to shift this distribution to higher redshifts, as predicted for example by the 
model of Rowan-Robinson (2001) shown. Chary et al (2004) give a predicted redshift distribution 
for the deeper (20 /tJy) GOODS-SV survey, based on a fit to the 24 /im counts, which has a broad 
peak at z = 0.8-2. 

Yan et al (2004) have discussed the nature of SPITZER 24 /im sources with reference a plot 
of S(24)/S(8) versus S(24)/S(r) and drawn attention to sources with high values of S(24)/S(r) as 
a potential new population. Specifically they find that 23 % of the FLS sample have S(24)/S(r) 
> 300 and suggest these are probably luminous infrared galaxies at z > 1. They find that 2% of 
their sample have S(24)/S(r) > 1000 and speculate that these may by obscured AGN at z > 0.6. 
In Fig 21 we show S(24)/S(r) versus redshift for the Lockman-VF with different symbols/colours 
for different sed types. For the objects plotted here, no new sed types are required. Galaxies with 
an infrared excess only at 24 /mi (open yellow triangles) and z > 1.5 are Type 1 AGN, according 
to their optical seds. Many of those at z < 1.5 could be Type 2 AGN. Figure 10 suggests that 
blank field sources (r > 25) are mainly AGN dust tori or M82 starbursts. A 24 //m optically blank 
source with S24 > 70 /iJy would have S(24)/S(r) > 200, so from Fig 21 would be likely to have 
z > 0.5. Whether the 24 /mi sources without optical counterparts represent a new population, or 
simply slightly more distant counterparts of the identified sources, remains to be seen. It is worth 
noting that S(24)/S(r) acts as a crude redshift indicator, in that if S(24)/S(r) < 5, then z < 0.5. 
However higher values can show a very wide range of z. 
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Fig. 20. — Redshift histogram for 24 fim sources brighter than 200 \xJy broken up by infrared sed 
type. Dotted curve: prediction of Rowan-Robinson (2001). 
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Fig. 21. — log 10 5(24)/5(r) for galaxies and quasars in Lockman-VF, with different symbols for 
infrared template types: filled circles (black): cirrus, filled triangles (red): M82 starburst, open 
squares (green): Arp220, crosses (blue): AGN dust tori, open circles (yellow): only one band with 
infrared excess. 
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8. Conclusions 

We have discussed optical associations, spectral energy distributions and photometric redshifts 
for SWIRE sources in the ELAIS-N1 area and the Lockman Validation Field. The band-merged 
IRAC (3.6, 4.5, 5.8 and 8.0 /xm) and MIPS (24, 70 and 160 /xm) data have been associated with 
optical UgriZ data from the INT Wide Field Survey in ELAIS-N1, and with our own optical Ugri 
data in Lockman- VF. Criteria for eliminating spurious infrared sources and for carrying out star- 
quasar-galaxy separation, and statistics of the identification rate, are given. 

The spectral energy distributions of selected SWIRE sources in ELAIS-N1 with spectroscopic 
redshifts, are modelled in terms of a simple set of galaxy and quasar templates in the optical and 
near infrared, and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst, 
and AGN dust torus) in the mid and far infrared. 

The optical data, together with the IRAC 3.6 and 4.5 data, have been used to determine 
photometric redshifts. For galaxies with known spectroscopic redshifts there is a notable improve- 
ment in the photometric redshift when the IRAC data are used, with a reduction in the rms scatter 
from 10% in (1+z) to 7%. While further spectroscopic data are needed to confirm this result, 
the prospect of determining good photometric redshifts for much of the SWIRE survey, expected 
to yield over 2 million extragalactic objects, is excellent. Some modifications to the optical-nir 
templates were required in the previously uninvestigated wavelength region 2-5 [im. 

The photometric redshifts are used to derive the 3.6 /xm redshift distribution. Agreement with 
model predictions are quite good. For those sources with a clear mid infrared excess, relative to 
the galaxy starlight model used for the optical and near infrared, the distribution of the different 
infrared sed types in the Li r /L opt versus Lj r plane, where Lj r and L opt are the infrared and optical 
bolometric luminosities, is discussed. Low luminosity cirrus galaxies have values of Li r /L opt similar 
to those found with IRAS data, but there is a new population of luminous cold cirrus galaxies with 
Lj r > L op t, which implies substantial dust extinction. These appear to present an enhanced rate of 
quiescent star formation compared to those seen in spirals in our neighbourhood. 

There is a surprisingly high incidence of galaxies with Arp220-like seds (10%), which was not 
predicted by the counts models of Rowan- Robinson (2001). They tend to have high values of 
Lir/L opt . There is also a high proportion of galaxies whose mid-ir seds are fitted by an AGN dust 
torus template (29%). Of these only 8 % of these are Type 1 AGN according to the optical-nir 
template fitting while 25% are fitted with galaxy templates in the optical-nir and have Lj r > L opt 
and so have to be Type 2 AGN. The remainder have Lj r < L op t so can be Seyferts, in which the 
optical AGN fails to be detected against the light of the host galaxy. The implied dust covering 
factor, > 75%, is much higher than that inferred for bright optically selected quasars. 

For galaxies detected at 70 or 160 /xm we have compared the infrared template fits and bolo- 
metric luminosity estimates derived from 4.5-160 /xm data, and from 4.5-24 /xm only. The rms value 
^<7io(Lf )0 ; j 4.5_i6o/-f'6oZ,4.5-24) is 0.30, corresponding to a factor 2 uncertainty in Lb Q i estimated from 
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4.5-24 (im only. This is unfortunate because we have 70 or 160 //m data for only 916 of the 126193 
SWIRE galaxies with photometric redshifts in Nl. Thus for 60 % of the SWIRE Nl galaxies with 
photometric redshifts, we have no information at all on the dust emission (no significant infrared 
excess relative to the starlight/QSO models), For 39 % we have an approximate estimate based 
on 4.5-24 /xm data which gives the bolometric luminosity only to an accuracy of a factor of 2, and 
for just 0.7 % of galaxies do we have a more accurate estimate based on longer wavelength data. 
Future instruments like PACS on Herschel, or large cooled apertures like SAFIR and SPICA, will 
be needed to fully quantify far infrared emission from galaxies at z > 1. 
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Table 1: Properties of SWIRE-EL AIS galaxies with plotted seds. The three sections of the table 
correspond to the galaxies of Figs 4-6, respectively, in order from the bottom of the figures. n se d 
is the infrared template type, and n typ is the optical template type (see text). Bracketed redshifts 
are photometric, rest are spectroscopic. 



SWIRE RA 


dec 


ELAIS name 


z 






C 




A v 


lcjL pt 




(2000) 








(Lo) 






(Lo) 




241.14116 


54.74209 


160433.8+544432 


0.0751 


1 


10.83 





2 





10.64 


241.47246 


54.37391 


150553.3+542225 


0.2116 


1 


11.23 


2 


2 





11.29 


241.90221 


53.95853 


160736.5+535731 


0.0298 


1 


10.08 


5 


3 





10.08 


241.90823 


54.76723 


160737.9+544601 


0.0913 


1 


10.94 


7 


2 





10.75 


242.01546 


54.88388 


160803.8+543202 


0.0528 


1 


10.83 


10 


2 





10.83 


242.05882 


54.47658 


160814.2+542836 


0.1198 


1 


10.87 


12 


3 





10.36 


242.14723 


53.83896 


160835.2+535022 


0.0627 


1 


10.48 


14 


1 





10.27 


242.21648 


54.21648 


160937.5+541258 


0.0862 


1 


10.74 


16 


1 





10.42 


242.76556 


54.72275 


161103.7+544322 


0.0627 


1 


10.33 


18 


3 





10.34 


243.28314 


54.86163 


161308.1+545141 


0.0700 


1 


10.38 


20 


3 





10.42 


241.67583 


55.83316 


- 


(0.479) 


1 


12.38 


-3 


4 





11.42 


242.62575 


55.02607 


161030.1+550135 


0.2740 


1 


11.78 





1 





11.19 


244.18948 


54.51995 


161645.3+543111 


0.2221 


1 


11.40 


2 


3 





11.80 


240.84489 


54.71045 


160322.8+544237 


0.2160 


2 


11.32 


5 


3 





10.78 


241.58780 


53.72672 


160621.0+534336 


0.2067 


2 


11.26 


7 


2 





11.13 


241.92168 


55.03159 


160741.1+550152 


0.4599 


2 


12.02 


9 


2 





10.66 


242.10521 


54.63612 


160825+543809 


0.9070 


2 


12.45 


11 


4 





11.40 


242.23936 


54.17452 


160857.6+541027 


0.2687 


2 


11.15 


14 


1 





10.88 


243.74185 


54.14808 


161457.9+540853 


0.2150 


2 


11.27 


16 


4 





10.93 


243.96436 


54.25998 


161551.3+541536 


0.2148 


2 


11.27 


18 


4 





10.93 


241.39261 


54.69358 


160534.3+544136 


(0.148) 


3 


10.90 


-5 


1 





10.44 


242.20854 


54.73938 


160850+544422 


(0.514) 


3 


11.96 


-3 


2 





11.32 


242.25496 


54.30231 


160901+541808 


0.3320 


3 


11.50 





4 


1.0 


10.80 


242.52415 


54.17478 


161005.8+541029 


0.0636 


3 


10.63 


3 


2 





10.84 


242.58041 


54.36514 


161019+542153 


0.2074 


3 


11.22 


6 


4 


0.7 


10.45 


241.73061 


53.67138 


160655.5+534016 


0.2136 


4 


11.13 


10 


7 





11.19 


242.30492 


53.90829 


160913+535429 


0.9924 


4 


12.14 


11 


7 





12.48 


243.15944 


53.38187 




2.138 


4 


12.98 


12 


7 





13.48 


241.77142 


53.59966 




3.653 


4 


13.80 


13 


7 





13.75 


240.75597 


54.75597 




0.7280 


4 


12.54 


17 


7 





12.29 



-38- 



Table 2: Statistics of SWIRE Catalogue in EN1 and Lockman-VF 



area 

total SWIRE Catalogue 
optically blank 
optical ID but no r, or not 4 bands 
(out of 3.6, 4.5 urn, UgriZ) 
stars 
phot z failed 
galaxies with phot z 

of which 4.5 //m 
5.8 fim 
8.0 /im 
24.0 /im 
70 /im 
160 



6.5 
301021" 
96155 6 
32721 

38500 
7452 
126193 

101139 

27554 
27074 
38601 
870 
406 



0.3 sq deg 
19140 a 
2991 
2331 

1154 
966 
11698 

8931 
2109 
1472 
878 
22 
14 



E 


29383 (23%) 


2361 (20%) 


Sab 


19326 (15%) 


1292 (11%) 


Sbc 


13483 (11%) 


1636 (14%) 


Scd 


29235 (23%) 3209 (28$) 




Sdm 


9512 (7%) 


949 (8%) 


sb 


19095 (15%) 


1769 (15%) 


QSO 


6159 (5%) 


481 (4%) 


cirrus 


8091 (31%) 


177 (25%) 


M82 


7595 (29%) 


166 (24%) 


A220 


2624 (10%) 


165 (24%) 


AGN dust torus 


7598 (29%) 


197 (28%) 


single band ir excess 


20793 


749 


no ir excess 


79492 


10244 



a all detected at 3.6 /im 
b includes small number of 
sources in WFS chip gaps 



